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Abstract 

A series of numerical simulations of Jet-A spray reacting flow in a single-element lean direct 
injection (LDI) combustor have been conducted by using the National Combustion Code (NCC). The 
simulations have been carried out using the time filtered Navier-Stokes (TFNS) approach ranging from 
the steady Reynolds-averaged Navier-Stokes (RANS), unsteady RANS (URANS), to the dynamic flow 
structure simulation (DFS). The sub-grid model employed for turbulent mixing and combustion includes 
the well-mixed model, the linear eddy mixing (LEM) model, and the filtered mass density function 
(FDF/PDF) model. The starting condition of the injected liquid spray is specified via empirical droplet 
size correlation, and a five-species single-step global reduced mechanism is employed for fuel chemistry. 
All the calculations use the same grid whose resolution is of the RANS type. Comparisons of results 
from various models are presented. 


1.0 Introduction 

A prerequisite for accurate prediction of turbulent combustion in the combustor is the ability of the 
turbulence model to capture the unsteady turbulent structures responsible for the mixing in the 
combustor. Approaches such as the large eddy simulation (LES) and the time filtered Navier-Stokes 
simulation (TFNS) are capable of capturing the dynamically important, unsteady turbulent flow 
structures. In the case of LES, the filtered equations are established by applying a spatial filter to the 
exact form of the governing equations. The filter width is typically the local grid size; in addition, the 
eddy viscosity has the local grid size as a model parameter. Therefore, the grid resolution and the model 
fidelity are formally linked, and, in principal, a grid independent solution cannot be cleanly reached. In 
the case of TFNS, the filtered equations are established by applying a temporal filter to the exact form of 
the governing equations. The filter width does not relate to the time step of the numerical solution, and 
the eddy viscosity contains the so called “resolution control parameter” which is conceptually defined as 
the ratio of the temporal filter width to a characteristic integral time scale. Since the grid resolution and 
the model fidelity are not formally linked, in principal, a grid independent solution can be unambi- 
guously attained. It should be pointed out that TFNS is not LES, nor hybrid RANS/LES, nor, in general, 
unsteady Reynolds-averaged Navier-Stokes (URANS). In our previous reports (e.g., Ref. 1), the TFNS 
was known as the partially resolved numerical simulation (PRNS) and the dynamic flow structure 
simulation (DFS) as the very large eddy simulation (VLES). In the present report, TFNS/DFS and 
PRNS/VLES are interchangeable. 
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Turbulent combustion also requires modeling the turbulence-chemistry interaction process. To this 
end, various sub-grid combustion models have been invoked, e.g., well-mixed model, eddy-break-up 
model (EBU) (Ref. 2), thickened flame model (Ref. 3), flamelet-based model (Ref. 4), conditional 
moment closure method (CMC) (Ref. 5), filtered mass density function/probability density function 
method (FDF/PDF) (Ref. 6), and linear eddy mixing model (LEM) (Ref. 7). In the well-mixed model, 
the distributions of scalars within a CFD computational cell are assumed to be homogeneous, i.e., their 
values equal to their respective mean values provided by the CFD solution of the filtered equations. The 
well-mixed model is often used in the direct numerical simulation (DNS). 

Fuel chemistry models ranging from global reduced mechanism involving a few species and reaction 
steps to detailed mechanism involving several hundreds of species and reaction steps exist. Global 
reactions are computationally less intensive, but they lack emissions and extinction information. Kinetics 
calculations using detailed mechanism consume long CPU time as well as large computer memory; 
hence, they are rarely, if ever, performed on the flight during the simulation. Typically, libraries of pre- 
computed tables are established by using detailed mechanism, then, linked to the simulation through 
interpolation procedure. Tabulation strategy such as the artificial neural network (ANN) has been used to 
further speed up the kinetics calculations while reducing the memory requirement (Ref. 8). 

For liquid fueled combustion, models for primary atomization, secondary droplet breakup, droplet 
vaporization, and droplet transportation are needed. The performance of these models directly impacts 
the local fuel vapor distribution in the calculated combustion field, therefore, can play a major role in the 
overall accuracy of the prediction. Last, but not the least, surrogate of the targeted fuel, e.g., the Jet-A, 
and its transport properties as well as its chemical kinetics should also be determined. 

The objective of this report is to summarize the calculated results from a series of numerical 
simulations of two-phase reacting flow in a single-element lean direct injection (LDI) combustor using 
the National Combustion Code (NCC). The NCC is a NASA in-house code for the modeling and 
simulation of multi -phase turbulent combustion in combustors (Ref. 9). The LDI concept has the 
potential for low emissions under operational (high temperature, high pressure) conditions. In this 
concept, the liquid fuel is injected from a venturi directly into the incoming swirling airflow, and the 
swirling airstream is used to atomize the injected liquid as well as to promote fuel-air mixing (Ref. 10). 
The flame structure can be very complex and locally range from non-premixed to premixed burning. 

This report is organized as follows. In the next section, a description of the TFNS formulation of 
two-phase flow and the candidate sub-grid combustion models are presented. In Section 3, some features 
of the experimental and computational setup are highlighted. This is followed by comparisons between 
measured data and computed results, and the concluding remarks. 

2.0 Formulation 

The conservation equations for compressible reacting flow are solved using the TFNS approach. To 
simulate spray combustion, Lagrangian droplet model is concurrently solved with the Eulerian gas flow. 

2.1 Gas Phase 

2.1.1 Governing Equations for Time Filtered Field 

Using the TFNS/DFS (previously known as PRNS/VLES) approach detailed in Reference 1, the 
following filtered equations are obtained, 

2. 1.1.1 Continuity and Momentum Equations 

P,f + (P“i), f =0 » P = P R f (!) 
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We call the ratio Aj/T 7 the resolution control parameter (RCP). It controls, conceptually at the governing 
equation level, the time scales content of the unsteady flow structures targeted for capturing in the 
practical numerical solution. 

When RCP ~ 1.0, all time scales of fluctuations (turbulence) have been filtered, and the TFNS 
simulation is a de facto RANS simulation, i.e., the directly calculated field is intrinsically the time mean, 
no turbulent fluctuation occurs in the directly calculated quantities. As the value of RCP decreases, the 
unsteady fluctuations will become more pronounced in the directly calculated field, if the grid sizes are 
adequate. To perform a TFNS simulation, we need to choose a value of RCP from the outset. The 
following relationship suggests a way to guide the selection of the value of RCP: 


RCP - 


k z 


^ ref / ^ref ^ ref 


( 6 ) 


For example, RCP = 0.3 means that we intend to directly calculate (i.e., to capture) flow structures 
associated with time scales responsible for about 70 percent of the total fluctuating (turbulent) kinetic 
energy, and the rest 30 percent of the total turbulent kinetic energy will be accounted for by the modeled 
sub-filter kinetic energy. 


2. 1.1.2 Transport Equations for Sub-Filter k-z 
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2. 1.1.3 Energy and Species Equations 


where 


dp e dp U i e _ dq t 


dt dxj 


- _ ^r x + PShk +P z + Q 


dx t 


dp® m _ + ?pUi® m 


dt 


dX; 


dx i 


<m) gpO;, 


dxi j 


+ P Sm Tfl — 1,2,***,7V 


P U^ m = pUi®m -r { T m) ^-(q 5,y + c 2 Qij ) 

p[/ z ^ has the same form as pt/^ , except that <P Z is replaced by e , and the 
diffusivity is replaced by turbulent heat conductivity. 
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2.1.2 Models for Sub-Grid Mixing and Combustion 

2. 1.2.1 Well-Mixed Model 

In well-mixed model, the distributions of scalars within a CFD computational cell are assumed to be 
homogeneous, i.e., their values equal to their respective mean values provided by the CFD solution of 
the filtered equations. Therefore, the filtered reaction rates are directly calculated using the filtered 
species mass fractions and filtered temperature defined at the computational cell center. The well-mixed 
model is, at the best, valid only in cells in which the sub-grid fluctuations are physically low. 

2. 1.2.2 Linear-Eddy Mixing (LEM) Model 

The LEM model (Ref. 7) is implemented in terms of a fractional splitting technique; it is divided into 
two processes: sub-grid and super-grid. The super-grid process emulates the convection of the scalar 
field by the grid-resolved velocity field across the surfaces of the computational cell. The sub-grid 
process, which occurs within each computational cell, consists of four operators: (a) molecular diffusion, 
(b) finite -rate kinetics, (c) volumetric expansion caused by the heat release, and (d) stochastic stirring 
due to the sub-grid eddies. 

2. 1.2.2. 1 Sub-Grid Process 

Within each CFD computational cell, a one-dimensional domain consisting of a fixed number of 
LEM elements is employed, the governing equations have the following form: 
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Species equation 


di i 
dt 


• m 

— -C+dYi-Yij+S^ 


Pe 


(13) 


Energy equation 

^L = -C^(h e -h) + S mle (14) 

F t stirring and F Y . st i rrmg represent the sub-grid turbulent mixing, and they are accounted for by 

employing a stochastic rearrangement of the LEM elements, known as the triplet mapping. This 
mapping, which mimics the characteristics of the stirring process inherent to turbulent flows, have three 
parameters: the stirring frequency X, a notional (stirring) eddy of size fi, and its stirring position on the 
linear LEM elements. The stirring position is randomly chosen from a uniform distribution. The eddy 
size is randomly selected from a prescribed distribution function given by 

/( / ) = (5/3)r 8/3 /[r| _5/3 -/r 5/3 ] (15) 

Here, the Kolmogorov length scale r| is estimated via r| ~l t Re -34 , and l, is the integral length scale of 
the turbulence. The event rate (frequency per unit length) is determined from 
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The stirring time interval between events is given by At siir = l/(AA) where A is the length of the 
one-dimensional LEM domain in a CFD computational cell. 

2.1.2.2.2 Super-Grid Process 

The grid-resolved convection of the scalar field is implemented by a Lagrangian transfer of LEM 
elements across the surfaces of the CFD computational cells. This Lagrangian transport is also known as 
splicing. Referring to Figure 1, for example, the outward mass through the right side of a cell computed 
from resolved velocity and density is equivalent to 1.5 LEM elements and colored in red. The outward 
mass through the bottom side of the mesh is equivalent to 2.5 LEM elements and colored in magenta. 
Similarly, the inward mass through the top and left sides of the mesh are equivalent to 6 LEM elements. 
Splicing will result in e.g., 14 LEM elements in this computational cell. In general, splicing will cause 
different computational cell to have different number of LEM elements. To avoid programming 
complexities in a parallel environment, the LEM domain is regridded to have the same fixed number of 
elements, and each element is of the same volume. Conservation of mass is maintained during 
regridding. 

In general, the LEM module receives the filtered velocity, pressure and turbulence fields from the 
CFD flow module. It also receives the source terms due to the liquid-phase contribution from the spray 
module. However, the feedback from the LEM module to the CFD flow module can be accomplished in 
several ways. In this report, the LEM module provides the species and the temperature fields to the CFD 
flow module and the spray module. Other feedback alternatives are currently being investigated. 
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2. 1.2.3 Filtered Mass Density Function (FDF/PDF) Model 

Under the present effort, the FDF/PDF model is implemented using an Eulerian Monte Carlo method 
(Refs. 1 1 and 12). The transport equation of FDF is solved by making use of an approximate 
factorization scheme. It is split into four operators respectively associated with the chemical reactions, 
the molecular mixing, the spray and the convection. In the Monte Carlo simulation, the filtered density 
function is represented by an ensemble of stochastic particles. Each particle carries enthalpy, temperature 
and species mass fractions. 

2. 1.2.3. 1 Sub-Grid Process 

For each notional particle, the following equations are solved: 

Species equation 

^XL = <-C^(Y i -Yi) + S mls (17) 

at 

Energy equation 

^L = -C^(h(-h)+S mle (18) 

where co = e/k , is an empirical constant. The first term on the right hand side of the energy equation 

accounts for the molecular mixing by making use of the relaxation to the ensemble mean. The last term, 
S m i e represents the contribution from the spray source terms. 

2.1.2.3.2 Super-Grid Process 

The grid-resolved convection of the scalar field is implemented in an Eulerian context. Referring to 
Figure 2, unlike the LEM splicing algorithm, the Eulerian convection and diffusion process is achieved 
through the content replacement of particles, as each Monte Carlo particle carries the information of the 
mass fraction of species, the enthalpy, the temperature, but not the volume. For example, the outward 
mass through the right side of a cell, computed from the grid-resolved velocity and density; and 
normalized by the mass in the cell, is equal to 5 percent and colored in red. The outward mass through 
the bottom side of the mesh is equal to 10 percent and colored in magenta. Similarly, the inward mass 
through the top and left sides of the mesh are equivalent to 40 percent and colored in dark and light 
green. If the total number of the particles is 20 in this cell, then one particle (5 percent of 20) will be 
randomly selected from its right adjacent cell, and its contents are copied to one randomly selected 
particle in the current cell. The same procedure is applied to other particles in colors. For those 
remaining particles (45 percent here) the contents are randomly shuffled within the cell. 

In general, the FDF/PDF module receives the filtered velocity, pressure and turbulence fields from 
the CFD flow module. It also receives the source terms due to the liquid-phase contribution from the 
spray module. However, the feedback from the FDF/PDF module to the CFD flow module can be 
accomplished in several ways. In this report, the FDF/PDF module provides the species and the 
temperature fields to the CFD flow module and the spray module. Other feedback alternatives are 
currently being investigated. 
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2.2 Liquid Phase 

The governing equations for the liquid phase are based on a Lagrangian formulation where the spray 
particle position and velocity are described by a set of ordinary differential equations. Various sub-models, 
such as the droplet drag model and the drop vaporization model, are needed to simulate the transport of a 
vaporizing spray particle. The specification of the fuel injector exit condition plays a major role in the 
fidelity of the simulation. Common practice is to specify the starting droplet condition using correlations of 
droplet sizes calibrated by relevant experimental data. In addition to the use of correlation, various models 
for primary atomization and secondary droplet breakup have been implemented into the NCC and 
evaluated/validated. A more detailed description of the liquid phase modeling and the coupling between the 
liquid and gas phase transport can be found in References 13 and 14. 

3.0 Experimental and Computational Setup 

The single-element LDI combustor is illustrated in Figure 3, more detailed description of the 
combustor geometry and the test rig can be found in Reference 15. Each element consists of an air 
passage with an upstream air swirler and a converging-diverging venturi section. The fuel is injected 
through the center of swirler and the fuel tip is at the throat of the venture. The air swirlers have six 
helical, axial vanes with downstream vane angles of 60°. The air then dumps into a combustion chamber 
with a square cross-section. Velocity measurements were taken with a two-component Laser Doppler 
Velociometry (LDV) system, temperature measurements were taken with thermocouples, and emissions 
data was gathered via an isokinetic probe and gas analyzer. Quartz makes up the combustion section. 

The combustor experiments may have significant convective and radiation heat losses. The temperature 
measurements reported are not corrected to adiabatic conditions. Experimental droplet measurements are 
collected with a Phase Doppler Particle Analyzer (PDPA). 

The computational domain consists of approximately 862,000 hexahedral elements (Fig. 3). At the 
air inflow boundary, the air flow speed is 20.14 m/s, the density is 1.19 Kg/m, and the static temperature 
is 294.28 K. At the combustor chamber exit, static pressure is specified for the steady RANS 
calculations, while convective conditions are used for the URANS and PRNS calculations. The operating 
pressure of the combustor is approximately 1 atm, and the measured pressure drop (as a percentage of 
the air inlet pressure) during the experiments was measured at 4 percent. 

The fuel is injected at 0.415 g/s, which gives a global equivalence ratio of 0.75. The specification of 
the starting condition for the fuel spray is particularly critical for accurate predictions. In this study, the 
following droplet size distribution is used (Ref. 13): 


— - 4.21 xlO 6 -±- 
n d i2 


-, 3.5 


l d 32 


0.4 


dd 

^32 


(19) 


Where n is the total number of the droplets and dn is the number of droplets in the size range between d 
and d+dd. This correlation also requires the specification of a Sauter mean diameter, <f 32 and the number 
of droplet classes. These specified inflow droplets will undergo evaporation without secondary breakup. 
Experimental data suggests a Sauter mean diameter around 32 pm, and the spray cone angle is 90°. 

In this study, the liquid fuel C^Ffy is used as the surrogate for the experimental Jet-A fuel, and a 
single-step, five-species global reduced mechanism (see e.g., Ref. 16) is employed for the chemical 
reactions Table 1. 
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4.0 Results 


The National Combustion Code (NCC) has been used for all the calculations. The results of the 
non-reacting flows and the comparison with the experimental data can be found in Reference 17, and 
they reveal the dominant flow structures in the LDI combustor. 

Figure 4 is a snapshot of the unsteady flow field. Embedded in this figure are the instantaneous 
iso-surface of the zero axial velocity component colored by the effective eddy viscosity and six 
instantaneous stream lines emanating from the upstream of the swirler, going through the converging- 
diverging nozzle, then passing through the combustion chamber. 

The dominant flow structures in the LDI combustor can be best visualized via the iso-surface of the 
zero axial velocity and the iso-surface of a relatively low pressure. The iso-surface of the zero axial 
velocity is also known as the vortex breakdown bubble (VBB). The iso-surface of a sufficiently low 
pressure captures the precessing vortex core (PVC). 

Figure 5 is a snap shot of the PVC and VBB. The dark blue region is a vortex core, which is formed 
near the venturi throat and extends into the combustor chamber. This spiraling vortex core rotates and 
breaks, it changes randomly in space and time. Embedded in this figure is an instantaneous stream line, 
which starts from the upstream of the swirler and goes through a complex, seemingly random path in the 
combustor chamber. This stream line spirals around the dark blue surface indicating that the dark blue 
region is indeed a vortex core. The light green surfaces are the iso-surfaces of the zero axial velocity. In 
addition to the VBB, there are some small structures near the dump plane and in the corner region. It is 
very reasonable to expect that the dynamics of the PVC and the VBB, as well as their interactions, are 
critical to the fuel-air mixing and the flame stability in the LDI combustors. 

The same grid is used for the two-phase reacting flow calculations. Figure 6 is a snapshot of the spray 
droplet field. In the following, we first report results from RANS, URANS, PRNS (all using the well- 
mixed model) and comparison with measured data, then the results from PRNS using different sub-grid 
models (LEM and FDF/PDF) and comparison with measured data. Due to the uncertainty in the imposed 
starting condition of the spray, the use of a simple single-step global chemistry model, and the yet to be 
optimized coupling between the CFD flow module and the LEM or FDF/PDF module in the framework of 
the time filtered Navier-Stokes simulation, the currently reported results and comparison with data should 
be viewed as preliminary. Results from refined simulations will be reported in the future. 

4.1 RANS, URANS, and PRNS Using Well-Mixed Model 

The steady state RANS (RCP = 1.0) solution is obtained by an iteration procedure using a Runge- 
Kutta scheme, and the static pressure is imposed at the outlet boundary of the computational domain. 
The unsteady URANS (RCP = 1.0) and PRNS (RCP = 0.26) solutions are obtained by a second order 
time accurate, iteratively implicit algorithm, and the outlet boundary condition facilitates the convection 
of pressure disturbances out of the computational domain. In this Section, the well-mixed model is 
employed to provide the filtered reaction source terms in the filtered species equations, i.e., the scalar 
field of the overall simulation is determined from the solution of the filtered energy and species 
equations of the finite- volume CFD module. 

Figure 7 presents the contours of the time averaged axial velocity in the center plane, while the time 
averaged temperature contours are shown in Figure 8. There are differences in the strength and extent of 
the center recirculation region and the corner vortex. More pronounced is the temperature difference in 
the flame zone. Figure 9 shows the comparison of the time averaged axial velocity along the center line. 
Comparisons of the time averaged axial velocity at several downstream cross-sections are presented in 
Figures 10 and 11. Comparisons of the time averaged temperature are given in Figures 12 to 14. 
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4.2 PRNS Using LEM and FDF/PDF Models 

For this set of calculations, the starting condition is an instantaneous solution of the PRNS using 
well-mixed model. Figure 15 depicts the starting condition in the center plane. With this starting field, 
the unsteady MacCormack scheme is chosen to continue the PRNS calculation using either LEM model 
or FDF/PDF model. When the LEM model is employed, 24 LEM elements per CFD computational cell 
are used, and up to 40 stirring events and reaction-diffusion operations are imposed over one CFD flow 
time step, rather than dynamically determined from the diffusion time scales and the correlation for the 
stirring events (see Section 2.1.2). When the FDF/PDF model is employed, 100 Monte Carlo particles 
per CFD computational cell are used, and a weighted time averaging scheme is adopted to reduce the 
statistical noise. In these two cases, unlike when the well-mixed model is employed, the scalar fields 
(i.e., temperature and species mass fractions) of the overall simulation are provided by the solution of the 
LEM module or the FDF module, while the velocity and pressure fields of the overall simulation are 
provided by the finite-volume CFD module. 

Figure 16 presents the contours of the time averaged axial velocity, temperature and fuel vapor mass 
fraction in the center plane obtained from LEM/PRNS. Figure 17 shows the time averaged axial velocity, 
temperature and fuel vapor mass fraction in the center plane obtained from FDF/PRNS. These two 
models result in quite different time averaged patterns. 

In the following, comparisons of time averaged results among LEM/PRNS, FDF/PRNS and the 
experiment are presented. The starting condition (denoted as well-mixed/PRNS) is also included for 
reference. The comparison of time averaged axial velocity along the center line is shown in Figure 18, 
the comparison of time averaged temperature along the center line is given in Figure 19. The time 
averaged temperature; axial velocity and azimuth velocity at several downstream stations are presented 
in Figures 20 to 22, respectively. Generally speaking, the comparisons do not suggest good agreement 
among the results. Some of the major contributing factors have been briefly mentioned above; the role of 
the CFD grid size in the performance of these sub-grid models also needs to be assessed. 

5.0 Concluding Remarks 

Numerical simulations of two-phase reacting flow in a single-element LDI combustor have been 
carried out using the time filtered Navier-Stotkes (TFNS) approach ranging from steady Reynolds averaged 
Navier-Stokes (RANS), unsteady RANS (URANS), to dynamic flow structure simulation (DFS). The 
TFNS/DFS is previously known as the PRNS/VLES. Results are obtained by employing different sub-grid 
mixing and combustion models, namely, well-mixed, LEM and FDF/PDF. All the calculations use the 
same grid whose resolution is of the RANS type. The present results should be considered as preliminary, 
because of the uncertainty in the imposed starting condition for the spray, the use of a five-species single- 
step global chemistry model, and the yet to be optimized coupling between the CFD finite-volume module 
and the LEM or FDF/PDF module. The role of the grid size in the performance of these sub-grid models 
also needs to be assessed. These issues are currently being addressed. 
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TABLE L— SINGLE STEP (GLOBAL) CHEMISTRY MODEL 



Reaction 

a 

(mole-cm-sec-k) 

n 

E 

(cal/mole) 

1 

4 C12H23+71 02=> 48 C02 + 46 H20 

GLO/C12H23 0.10/ 

GLO/02 1.65/ 

8.60E+1 1 

0.00 

3.00E+4 
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Arrows -- Convection direction 
Initial number of lumps - 12 
Light Green - inward mass (1 .5) 

Dark Green - inward mass (4.5) 

Red - outward mass (1 .5) 

Magenta - outward mass (2.5) 

Current number of lumps -- 14 

Final number of lumps will became 12 after re-gridding 
via equalizing the volume. 


Figure 1 . — Schematic illustrating the splicing algorithm used for scalar convection. (Assuming number of 
LEM elements is 12. Each element carries its own volume and species density). 



Arrows -- Convection direction 
Light Green - inward mass (15%=3) 

Dark Green - inward mass (25%=5) 

Red - outward mass (5%=1 ) 

Magenta - outward mass (1 0%=2) 

Number of Monte Carlo particles - 20 
Contents of randomly selected particles are 
replaced by those of randomly selected particles 
from adjacent cells. 


Figure 2. — Schematic illustrating the process of the scalar convection and diffusion for the Eulerian Monte Carlo 
FDF/PDF model. 


A Swirl-Stabilized LDI Flame Tube 



Grid Distribution for the LDI 


" • 001 OCT 

Combustor (861 823 hexahedron elements) 


Figure 3. — Swirler geometry and computational domain for the LDI combustor. 
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Figure 4. — A snapshot of the flow field from the PRNS simulation. 



Figure 5. — A snapshot of the flow structures: precessing vortex 
core (PVC) and vortex breakdown bubble (VBB). 
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Very large eddy simulation PRNS/VLES of the Jet-A spray 
turbulent combustion in a single element LDI combustor 



Figure 6. — A snapshot of the spray droplet field. 
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RANS, spray & reacting, cfl=1.0, chem_solver=0 
no_stream=8, no_group=20, flowf=4. 1 5E-4 
dtml=1.E-06, dtgl=5.E-04, dtil=5.E-04 
460,000 iterations 
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Mean of URANS, spray-reacting, int_chem_solver=-1 
no_stream=8, no_group=20, flowf=4.15E-04 
dtml=2.E-07, dtgl=2.E-06, dtil=2.E-06 
dt=2.0E-06, 2nd=-0.02, 4th=0.06, cfl=1.0 
73,000 time step (av over last 10,000 steps) 
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Mean of PR NS, spray-reacting, cfl=1 .0 
no_stream=8, no_group=20, flowf=4.15E-04 
dtml=1 E-07, dtgl=1 .E-06, dtil=1 .E-06, dt=1 .0E-06 
2nd=0.0, 4th=0.01 
1 10,000 time step 
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59 

49 

39 

29 

18 
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■I - 2 

H ' 2 
™ -22 



Figure 7. — Contours of the time averaged axial velocity in the 
center plane (well-mixed model): RANS, URANS, and 
PRNS. 
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RANS, spray & reacting, cfl=1 .0, chem_solver=0 
no_stream=8, no_group=20, flowf=4.1 5E-4 
dtml=1.E-06, dtgl=5.E-04, dtil=5.E-04 

460,000 iterations 


Y 



x 
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T 

2278 

2113 
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1615 

1450 

1284 

1118 

952 

787 

621 

455 



Mean of URANS, spray-reacting, int_chem_solver=-1 
no_stream=8, no_group=20, flowf=4.15E-04 
dtml=2.E-07, dtgl=2.E-06, dtil=2.E-06 
dt=2.0E-06, 2nd=-0.02, 4th=0.06, cfl= 1 .0 

73,000 time step (av over last 10,000 steps) 
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Mean of PRNS, spray-reacting, cfl=1 .0 
no_stream=8, no_group=20, flowf=4.15E-04 
dtml=1 .E-07, dtgl=1 .E-06, dtil=1 .E-06, dt=1 .0E-06 
2nd=0.0, 4th=0.01 

1 10,000 time step 
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1300 
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Figure 8. — Contours of the averaged temperature in the 
center plane (well-mixed model): RANS, URANS, and 
PRNS. 
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Figure 9. — Comparison of the time averaged axial velocity along the center line (well-mixed model): RANS, URANS, 
and PRNS. 
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Figure 10. — Comparison of the time averaged axial velocity at x = 9 mm (well-mixed model): RANS, URANS, 
and PRNS. 
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Figure 1 1 . — Comparison of the time averaged axial velocity at x = 29 mm (well-mixed model): RANS, URANS, and 
PRNS. 
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Figure 12. — Comparison of the time averaged temperature along the center line (well-mixed model): RANS, URANS, 
and PRNS. 
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Figure 13. — Comparison of the time averaged temperature at x = 5 mm (well-mixed model): RANS, URANS, 
and PRNS. 
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Figure 14. — Comparison of the time averaged temperature at x = 10 mm (well-mixed model): RANS, URANS, 
and PRNS. 
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Axial velocity, m/s 
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C12H23, mass fraction 


0 0.1071 0.2143 0.3214 0.4285 0.5356 



Figure 15. — Starting condition for LEM/PRNS and FDF/PRNS (in the center plane). 
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Temperature, K 
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C12H23, mass fraction 

0 0.1071 0.2142 0. 3214 0.4285 0.5356 


Figure 16. — LEM/PRNS time averaged contours in the center plane: axial velocity, temperature, and 
fuel vapor mass fraction. 
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Figure 17. — FDF/PRNS time averaged contours in the center plane: axial velocity, temperature, and 
fuel vapor mass fraction. 
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Figure 18. — Comparison of time averaged axial velocity along the center line 
(the instantaneous starting condition denoted as well-mixed/PRNS is included 
for reference). 



Figure 19. — Comparison of time averaged temperature along the center line 
(the instantaneous starting condition denoted as well-mixed/PRNS is included 
for reference). 
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Figure 20. — Comparison of time averaged temperature at x = 5 mm, 10 mm 
(the instantaneous starting condition denoted as well-mixed/PRNS is 
included for reference). 
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Figure 21 . — Comparison of time averaged axial velocity at x = 9 mm, 20 mm 
(the instantaneous starting condition denoted as well-mixed/PRNS is 
included for reference). 
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Figure 22. — Comparison of time averaged azimuth velocity at x = 9 mm, 20 mm 
(the instantaneous starting condition denoted as well-mixed/PRNS is included 
for reference). 
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